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The formation of a higher-order singularity by the coalescence of two exceptional points 
is studied in the anti-hermitian limit of the complex-extended 3-level Richardson-Gaudin 
model 



1. Introduction 

An essential part of the motion of short-lived nucleonic matter is in classically 
forbidden regions and, henc e, its properties are effected by both the continuum 
and many-body correlations uJ&l. The effect of resonances and the non-resonant 
scattering states can be considered in the OQS extension of the shell model (SM), 
the so-called continuum shell model (CSM) H Two realizations of the CSM have 
been studied recently: the real-energy CSM ^ and the complex-energy CSM EEl j 
the so-called Gamow Shell Model (GSM). For hermitian Hamilton operators, both 
real-energy CSM and complex-energy CSM (GSM) lead to the complex-symmetric, 
non-hermitian eigenvalue problem. As a result, OQSs exhibit unintuitive properties, 
qualitatively different from closed quantum systems (CQSs). 

One of them is the phenomenon of sharp crossing of resonances with same 
quantum numbers (same symmetries). Among degeneracies associated with avoided 
crossings in quantal spectra, one finds a diabolic point (DP) El, and an exceptional 
point (EP)^, which appears in complex extended hermitian Hamiltonians. In this 
Letter, we shall discuss some aspects of degeneracies in spectra of a prototypical 
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OQS, the 3-level Richardson-Gaudin (RG) model. 

2. The model 

RG models IS are based on the SU(2) algebra with elements Kp, K^, and K®, 
fulfilling the commutation relations: [Kf, Kf] — 8wK^ , [Kf, Kp] = ±6wKp , 
where indices I, V refer to a particular copy from a set of L, SU{2) algebras. In the 
following, we shall use the pair representation of SU(2) algebra. The elementary 
operators in this representation are the number operators Nj and the pair operators 
A\, j4», defined as: 

N S = £ a] m a jm ; A) = £ a] m a) m = (A^ (1) 

TCI TCI 

where j is the total angular momentum and m is the z-projection. The state jfn 
is the time reversal of jm. The relation between the operators of the pair algebra 
and the generators of the SU(2) algebra is: 

K? = ijyr I -in I ; k+ = (k^ = \a\ (2) 

where is the particle degeneracy of level I. With this correspondence, one can 
introduce the non-integrable 3-level pairing Hamiltonian: 

H(g)=J2^Ni+gJ2 A i A j+9'J2 N i ' 9' = 19 (3) 

i ij i 

where 7 < is a C- number and coupling constants g, g' are complex. For 7 = 0, 
the Hamiltonian §5§ is integrable. 

2.1. Level degeneracies 

The position of degeneracies in the complex g-plane are indicated by the roots of 
coupled equations: 

det[H(g)-EL]=0; ^det [H (g) - EI] = (4) 

Eliminating E from these two equations, one is left with a polynomial in g of degree 
M < n(n — 1). The maximal degree of D(g) (M = n(n — 1)) appears for non- 
integrable Hamiltonians. The eigenvalue degeneracies E m (g) — E m i(g) at g = g a 
[a. = 1, . . . , M), can be found numerically by looking for sharp minima of D(g) in 
the complex (/-plane. Degenerate eigenvalues can be single-root (EP) or double-root 
solutions of D(g). 

3. Formation of the pseudo-diabolic point in the anti-hermitian 
Hamiltonian 

In the following, we shall analyze solutions of the 3-level non-integrable pairing 
model ([3]) in the case of 2 pairs of fermions in a valence space of degeneracy £li = 
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Fig. 1. The formation and decay of a pseudo-DP at g' = —g/2 (the blue star) in the anti-hermitian 
limit of the non-integrable 3-level pairing Hamiltonian j5"|l. 
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Fig. 2. Real part of phases for all eigenvectors along the trajectory surrounding the pseudo-DP 
at g = — i/(4v2). The phase of each vector is defined with respect of its reference value at <j> = 0. 
For more details, see the discussion in the text. 



2, 17 2 = 6, ri 3 = 2. Energies of levels are E\ — 0, £2 = 1, and £3 = 2. In this model 
space, there are n = 4 many-body states, and the total number of roots of D(g) 
is M = 12. The 3-level pairing model is reflection symmetric with respect to the 
real g-axis (Tm(g) = 0). In the lower half-plane, all eigenvalues are either discrete 
states at TZe(g) = or decaying resonances. Complex conjugate states are situated 
in the upper half-plane (Tm{g) > 0) and correspond to capturing resonances. 

Fig. [1] shows a pattern associated with the formation/breakup of the pseudo- 
DP at g = — i/(4\/2), i.e. in the non-hermitian limit of the model ([3|). With an 
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Fig. 3. Real and imaginary parts of eigenvalues at around the pseudo-DP g = — i/(4%/2) are plotted 
along the cut (0,Zm(g)) in (/-plane. For more details, see the description in the text. 

increasing value of g' /g, the two EPs in the complex g-plane approach each other, 
coalesce at g'/g = —1/2 and move along the axis IZe(g) = for g'/g > —1/2. At 
g = — i / (4V2) for g'/g = —1/2, real parts of all eigenenergies are identical. 

Fig.[2]shows a real part of the phase along the trajectory in g-plane surrounding 
the pseudo-DP, for all eigenvectors in the chosen model space. The phase of each 
vector is defined with respect of its value at <j> — 0. The trajectory is defined by: 
g(<j>) = -(l/(4V2)i + O.Ole'f The (complex) eigenvalues for vectors '2' and '3' 
are degenerate at the pseudo-DP and change their phase by 7r at each encircling, 
as expected for an ordinary DP. However, in contrast to the ordinary DP, both 
eigenvalues and eigenvectors merge forming a Jordan block, like for an ordinary 
EP. Nevertheless, the crossing of eigenvalues is sharp as in the case of ordinary 
level crossing or the DP. One can see in Fig. [2] that phases of vectors T' and '4', 
which have different imaginary parts of eigenvalues, remain approximately constant. 

The complex behavior of eigenvalues in the neighborhood of a pseudo-DP in 
complex g-plane are shown in Fig. [3] along the cut (0,lm(g)). In the lower part, 
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one can see the appearance of a pseudo-DP at g = -i/(4V2) ~ -0.177i as a 
crossing of imaginary parts of eigenvalues '2' and '3'. At g ~ — i/(4v2) 5 real parts 
of all eigenvalues are equal (see the upper part of Fig. [3]). Notice at g = — 0.142276i, 
close to the pseudo-DP, a formation of an EP by eigenvalues '1' and '3'. 

4. Outlook 

In contrast to a CQS (a hermitian system), the exact degeneracies of eigenvalues 
for states of the same quantum numbers (symmetries) are allowed in an OQS. Such 
degeneracies (EP, D P, pseudo-DP) are closely related to topological anomalies of the 
Hilbert space "^^l j n this Letter, we have demonstrated the formation of a pseudo- 
DP in the anti-hermitian limit of a 3-level pairing Hamiltonian. It remains to be 
seen in future studies using the CSM and realistic two-body interactions whether 
those anomalies can have observable consequences in resonance spectroscopy and 
low-energy nuclear reactions. 
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